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The Shockley equation (SE), originally derived to describe a p-n junction, was frequently used in the past 
to simulate current-voltage (j/V) characteristics of organic solar cells (OSC). In order to gain a more de- 
tailed understanding of recombination losses, we determined the SE parameters, i.e. the ideality factor and 
the dark saturation current, from temperature dependent static jA'-measurements on poly(3-hexylthiophene- 
2,5-diyl)(P3HT):[6,6]-phenyl-C6i butyric acid methyl ester (PCBM) bulk heterojunction solar cells. As we 
show here, these parameters are directly related to charge carrier recombination and become also accessible 
by transient photovoltage and photocurrent methods in the case of field-independent charge carrier generation. 
Although determined in very different ways, both SE parameters were found to be identical. The good agree- 
ment of static and transient approaches over a wide temperature range demonstrates the validity of the Shockley 
model for OSC based on material systems satisfying the requirement of field-independent polaron-pair disso- 
ciation. In particular, we were able to reproduce the photocurrent at various light intensities and temperatures 
from the respective nongeminate recombination rates. Furthermore, the temperature dependence of the dark 
saturation current jq allowed determining the effective band gap of the photoactive blend perfectly agreeing 
with the literature values of the energy onset of the photocurrent due to charge transfer absorption. We also 
present a consistent model directly relating the ideality factor to recombination of free with trapped charge car- 
riers in an exponential density of tail states. We verify this finding by data from thermally stimulated current 
measurements. 
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I. INTRODUCTION 

Organic solar cells (OSC) hold the potential to become a 
low-cost photovoltaic technology, produced in a roll-to-roll 
process and thus deserve serious consideration.^ Therefore, 
the efficiency of organic solar cells still has to be improved 
and a proper understanding of the elementary physical pro- 
cesses is required. Charge carrier recombination is one of the 
particularly crucial processes as it limits the electrical perfor- 
mance of OSC. Indeed, the polaron recombination dynamics 
fed into the continuity equation have been shown to deliver 
sufficient information to reproduce the complete j/V response 
of P3HT:PCBM solar cells at room temperature.^ Recently, 
Maurano et al^ successfully predicted the open circuit volt- 
age of OSC based on four different polymer: fullerene compo- 
sitions by studying nongeminate recombination. Thereby, au- 
thors combined parameters obtained by static (jA') and tran- 
sient photovoltage / transient photocurrent (TPV/TPC) meth- 
ods in the Shockley equation.'* In order to justify the valid- 
ity of this approach and/or its limitations, we i) directly com- 
pared the ideality factor and dark saturation current density 
derived by the static and transient measurements and ii) re- 
produced the photocurrent for P3HT:PCBM OSC over the 
temperature range from 200 to 300 K at various light inten- 
sities. We demonstrate that identical Shockley parameters can 
indeed be determined from transient and static measurements 
making SE applicable to OSC, if the charge carrier generation 
is voltage independent. Deviations found at lower tempera- 
tures are attributed to voltage dependent charge canier photo- 
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generation due to the lack of thermal activation energy. Ex- 
tending the work of Kirchartz etal.,^ we analyzed the impact 
of charge carrier recombination on the ideality factors in de- 
tail, and derived the dominant nongeminate loss mechanism. 

n. THEORETICAL BACKGROUND 

The established model to describe the current-voltage re- 
sponse of semiconductor p-n junctions under illumination is 
based on the ideal SE, 

with dark saturation current jo, elementary charge q, ideal- 
ity factor «,v/, thermal energy kT and the photogenerated cur- 
rent jgen-'^ Conventional p-n junctions are characterized by 
well developed energy-bands where photogeneration yields 
delocalized charge carriers. In contrast, OSC are based on 
excitonic materials, in which the photogeneration of quasi 
free polarons upon photon absorption is via bound precur- 
sor states and strongly depends on the active material. Fur- 
thermore, BHJ devices consist of donor and acceptor material 
phases across the whole volume, resulting in many spatially 
distributed he teroj unctions instead of a single planar one. 

Despite these differences to inorganic p-n junctions, the 
SE was successfully applied to organic solar cells based on 
different donor-acceptor systems in the past.^^ Their appli- 
cation to devices based on MDMO-PPV: PCBlVpl failed due 
to voltage dependent photocurrent generation jge„{V). For 
OSC based on P3HT:PCBM, used in this work, the polaron 
pair dissociation is reported to be independenPHSl or weakly 
dependen^ii^s^on voltage, providing a good starting point for 
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our analysis. For simplicity we assumed an ideal diode in this 
study, neglecting any series resistance (R^ = 0) or leakages 
by shunts (Rp = °o). No significant influence on the data was 
found by taking a series resistance into consideration, which is 
experimentally determined from the ohmic range of dark j/V- 
characteristics. Eq. ([TJ solved for Voc at a given light intensity 
with j{V„c) = and for jg^,, 3> jo results in 

K.=n,./— Inf^ + l) (2) 
1 \ .10 ) 

«n,y — In^- (3) 
q JO 

We note that in order to experimentally access the photogen- 
erated charges by jA' measurements, we assumed a voltage 
independent generation current jge,,. 

III. EXPEMMENT 

Sample preparation 

Bulk heterojunction (BHJ) solar cells were prepared 
by spin coating a 35 nm thick layer of poly(3,4- 
ethylendioxythiophene):polystyrolsulfonate (Clevios P VP AI 
4083) on indium tin oxide samples with post-annealing step 
of 130°C for 10 minutes in a water-free environment. The 
P3HT:PCBM 1:0.8 blend made from solutions of 25 mg/ml 
in chlorobenzene was spin coated in an inert atmosphere re- 
sulting in 100-120 nm thick active layers. After annealing 
for 10 min at 130°C Ca (3 nm)/Al (100 nm) were evaporated 
thermally on top of the organic layer. P3HT was purchased 
from Rieke Metals (P200, > 98 % regioregular), PCBM (pu- 
rity 99.5 %) from Solenne. All materials were used without 
further purification. Prior to any additional measurements an 
Oriel 1 160 AM1.5G solar simulator was used to perform illu- 
minated j A'-measurements of devices kept under inert glove- 
box atmosphere. The BHJ solar cells showed a power con- 
version efficiency of approximately 3.3 %. The samples were 
transferred into a closed cycle optical cryostat for temperature 
dependent static and transient electrical studies. 

Experimental techniques 
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FIG. 1. (Color Online) (a) Charge carrier density and (b) small 
perturbation charge carrier lifetime vs. open circuit voltage from 
TPV/TPC measurements for a solar cell based on P3HT:PCBM in 
dependence on temperature, as indicated by arrows. Voc was varied 
by changing the illumination intensity. 



1.5 Gi2 impedance adapter to ensure open circuit conditions. 
Both, voltage and current transients were acquired by a digital 
storage oscilloscope (Agilent Infinium DSO90254A). 



IV. RESULTS AND DISCUSSION 

The results of temperature dependent TPV/TPC measure- 
ments are depicted in Fig. [T] As can be seen, charge car- 
rier density n (Fig.[T](a)) and small perturbation lifetime Tah 
(Fig.[T](b)) both exponentially depend on Voc and can be de- 
scribed by 



TPV measurements monitor the photovoltage decay upon a 
small optical perturbation during various constant bias light 
conditions. The voltage transient enable to determine the 
small perturbation charge carrier lifetime Tajj in dependence 
of the respective open circuit voltage due to the constant back- 
ground illumination. TPC is measured with identical optical 
perturbation, but under short circuit conditions. Further ex- 
perimental details can be found in Ref. [H [U TPV/TPC 
measurements were performed using a 10 W white light LED 
(Seoul) and focusing optics for bias light illumination and an 
attenuated Nd:YAG laser flash (A- = 532 nm, < 80 ps pulse 
duration) to generate an additional small amount of charges 
within the device. The voltage transient was probed using a 



with n„ defined as the ideality factor of charge carrier density. 
And similarly, 

TA„=TA„oexp(^-— j (5) 

with the ideality factor of charge carrier lifetime n-^. We note 
that the slopes in the semi-logarithmic representation of Fig. [T] 
were previously described by parameters y and pjSl which are 
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(6) 
(7) 



related to our dimensionless ideality factors n„ and «x by 

- g 

q 

As seen in Fig.[T]both slopes (P,Y) increase when the temper- 
ature is lowered. 



Equivalence of SE parameters from static and transient 
metliods 



To describe the experimentally found polaron dynamics in 
organic photovoltaic devices a generalized equation dn/dt = 
-kxn^+^ = -R{n) is often used.^^^From TPV/TPC analy- 
sis, the decay order of the recombination rate R,'k+\, and the 
small perturbation charge carrier lifetime lAn can be experi- 
mentally determined. Using these values the total charge car- 
rier lifetime i„ can be calculated, as it was shown in Ref.lTTl 



^71 = 'i:An(A.+ 1), 



(8) 



Using Eqs. (jsj) and (|8j), small perturbation life- 
times values can be eliminated. 



where A- — 



- T„o exp 



n^kT 

Dividing Eq. Q by Eq. (|9]l, we obtain: 



n 
1„ 
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where R{n) is defined as the recombination rate under illumi- 
nation, depending on illumination intensity and /?o the recom- 
bination rate in dark. Solving Eq. ( [T0| for V^c yields a more 
generalized expression. 
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with the recombination ideality factor tiR defined as 

-1 



(11) 



(12) 



which is consistent with earlier representations of the ideality 
factor.^ 

Comparing Eq. ( [TT] ), derived by using relations empirically 
found in transient measurements (Eqs. (|4]i,(|5]l), and the right 
hand side of Eq. the SE solved for Vac to describe steady- 
state j/V-measurements, the same equation structure becomes 
apparent. To compare both equations, first the equality of re- 
combination rates and respective currents is motivated. At 
Voc generation and recombination rates are equal, i.e. G = R, 
which implies the generation current (see Eq. ([3]l) is can- 
celled by the recombination current defined as 
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FIG. 2. (Color Online) (a) Ideality factor of a P3HT:PCBM 1:0.8 so- 
lar cell in dependence on temperature revealed from TPV/TPC stud- 
ies on the recombination dynamics and from illuminated j A' data 
analyzed with the respective Shockley equation (see Eq. |3](). (b) 
Respective dark saturation current densities in dependence on tem- 
perature, jiosso calculated from TPV/TPC studies on the recom- 
bination dynamics in the dark with Eq. \14\ and 70 from illuminated 
j/V response at Voc via Eq. l[3j. 



with the thickness of the active layer d. Likewise, the dark 
saturation current jo may be treated as thermally generated 
intrinsic charge carriers with the density «o and lifetime Tawq 
canceled by recombination in thermal equilibrium. Hence, the 
respective loss current in the dark y'/o„„ can be written as 



Jlo 



: Roqd. 



(14) 



jioss =R{n)qd, 



(13) 



The above considerations provide a rationale to the in- 
terconnection between the SE parameters derived diversely. 
Fig.[2](a) shows ideality factors n,v/ ™d hr, as defined above, 
which were experimentally obtained as function of tempera- 
ture by means of static (j/V) and transient (TPV/TPC) meth- 
ods, respectively. To determine the static SE parameters j/V 
characteristics were recorded for different bias lights. Then, 
Voc was plotted vs. ige„ and analyzed by Eq. Q to yield the 
ideality factor n,-^. Within the measured temperature range 
and «,v/ correspond well to each other. Similary, the dark sat- 
uration current jo and the dark recombination current jiossQ, 
depicted in Fig.|2](b), are almost equal for different tempera- 
tures. The ideality factor rises from 1 .2 at room temperature to 
1.5 at 170 K. A similar trend with temperature was observed 
in the past for organic solar cells based on different mater ials, 
although a stronger temperature dependence was obtained.'^'^ 

From the experimental observation shown in Fig. 2] the SE 
equivalence for P3HT:PCBM based OSC to Eq. (fHlTderived 
from studying charge carrier recombination by transient tech- 
niques, can be suggested. 



FIG. 3. (Color Online) Dark saturation current evaluated by study- 
ing the charge carrier dynamics plotted versus 1 /Tn{T) in order to 
estimate the effective band gap energy Eg from Eq. ijTSji. 



Effective band gap energy Eg 

We first focus on the strong temperature dependence of the 
dark saturation current (Fig. |2](b)). In Fig.|3] the dark satu- 
ration current jiosso determined by TPV/TPC is plotted versus 
I / {TnR{T)), with the temperature dependent ideality factor 
nR{T) from Fig.[2](a). The linear fit of the semi-logarithmic 
plot is in excellent agreement with the measured data and is 
expressed by the relation 
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(15) 



The fit demonstrates that the dark saturation current follows 
the Boltzmann distribution and allows to identify the effec- 
tive band gap energy Eg, which is proportional to the dif- 
ference between HOMO and LUMO energy levels of donor 
and acceptor, respectively. We found Eg — 1.07 eV for 
P3HT:PCBM. From the Voc{T) data (not shown), a value of 



0.9 eV was evaluated according to Ref. 



Vandewal 



et. al obtained the value of Eg — 1 .08 e V from the energy on- 
set of the photocurrent generated by CT absorption,^' which 
agrees particularly well with our value. Guan et al. found 
somewhat larger values {Eg « 1.3—1.4 eV) using ultravio- 
let and inverse photoemission spectroscopy.l^^As discussed in 
Ref.|23] energy values derived by photoemission spectroscopy 
can be seen as onsets of respective valence and conduction 
band. Due to charge carrier relaxation into the tail states of 
the density of states, the techniques such as TPV/TPC and j/V- 
measurements are expected to reveal the effective bandgap Eg 
from the relaxed density of occupied states instead. 



Loss current vs. photocurrent 
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FIG. 4. (Color Online) Current density jsat (solid lines) from j/V re- 
sponse compared to recombination current jioss (dashed lines) mea- 
sured by TPV/TPC for different temperatures and bias lights accord- 
ing to Eq. l |16[ ). For T <200 K a stronger deviation of both parame- 
ters becomes apparent, indicated by the dashed vertical line. 



If the generation of charge carriers is voltage independent 
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and 



for P3HT:PCBM cells in the range of - 1 V < V < K 
jam is not reduced by nongeminate recombination, we expect 



]mt — J I OSS- 



(16) 



For comparison, Eq. ( pjj ) was used to calculate the recombi- 
nation current //„ss under illumination at open circuit. 

The saturation current and the nongeminate recombi- 
nation current y'/oii at Vac determined via TPV/TPC measure- 
ments for different illumination intensities from 0.1 sun to 
about 1.8 suns are shown in Fig.|4] Within experimental er- 
ror, jioss is in a very good agreement with jsat for all light 
intensities and temperatures above 200 K. Below this temper- 
ature, jsat becomes smaller than y'/„ss, which seems to imply 
more charges recombining than being generated. This coun- 
terintuitive result can at least partly be explained by estimating 
the nongeminate recombination losses at V = — 1 V according 
to Ref. 24 We find 2% loss at room temperature and 5% at 
170 K, revealing the limitations of the analysis as the gener- 
ation current is underestimated. Accounting for the nongem- 
inate losses, a very small deviation of and jioss remains, 
showing that the field dependence of the photogeneration is 
less than 10% between open and short circuit conditions even 
at low temperatures. In accordance with earlier findings,^^^^ 
the photocurrent in P3HT:PCBM solar cells is at best slightly 
field dependent in this voltage range. 



To complete the comparison of Eq. ([3]l and Eq. ( 1 1 



we con- 
sider the photogenerated charge carrier density, being equal 
to the density of charges recombining at Voc- In order to es- 
timate the photogenerated charge carriers, we used the satu- 
rated current density under illumination in reverse bias 

(y = -iv). 



Contributions to recombination ideality factor 

In Eq. ( [TT| the recombination ideality factor hr was de- 
fined. As mentioned above riR is composed by the ideality fac- 
tor for charge carrier density n„ and charge carrier lifetime n^. 
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FIG. 5. (Color Online) (a) Temperature dependence of the inverse 
ideality factors for charge carrier density, n^^ , and lifetime, n^', as 
well as their sum, resulting in the inverse ideality factor (b) 
Corresponding recombination parameter 'k = nn/n^ (left) and recom- 
bination order A. + 1 (right). 



respectively (see Eq. ( 12 1). In Fig.|5](a), the inverse ideality 



factors derived from the transient experiments are plotted ver- 
sus temperature. At 300 K, we found «i: = 1 .84 and n„ — 3.49, 
yielding a recombination ideality factor of tin ~ 1-20. As 
shown already in Fig. |2](a), the ideality factor from the j/V- 
characteristics under illumination was nui = 1.21. Although 
we focussed on the ideality factors measured under illumina- 
tion, we state for comparison that the dark ideality factor was 
1.30, whereas Kirchartz etal.^ found a value of 1.6. 

The recombination order, given byA-+l = ^ + las pointed 
out above, is shown in Fig.|5](b). The recombination order at 
300 K of A, + 1 = 2.90 was slightly higher as compared to the 
values of 2.6 (Ref.QiJ or 2.75 (Ref.i5). 



Interpretation of ideality factors: Recombination mechanism 

Based on the very concise approach offered by Kirchartz 
etal.l^ we used the experimentally determined ideality fac- 
tors to consider the dominant recombination mechanism in the 
P3HT:PCBM solar cells under open circuit conditions. 

Ideality factors higher than 1 in organic bulk heterojunc- 
tion solar cells are commonly seen as an evidence for a dom- 
inant trap assisted recombination process. Kirchartz etal. 
proposed the recombination of free charge carriers tic with 
trapped charge carriers n, in an exponential tail as dominant 
nongeminate loss mechanism. The corresponding recombina- 
tion rate is R<>^ ricHt, assuming symmetric electron and hole 
concentrations. Using the relation, oc rtj of free n^- to 
trapped n, charge carriers they were able to calculate the char- 
acteristic (Urbach) tail energy Eu from experimental datP 
by three routes: (i) By the voltage dependence of the ex- 
tracted charge carrier concentration, which — rewritten in our 



notation — is given by the ideality factor n„, 

n„kT 



(17) 



using their experimental n„ w 4.2. (ii) From the (in their 
case dark) ideality factor, based on the relation by van Berkel 
etal.^ 



'id 



kT 
2Eu 



(18) 



which corresponds to our Eq. ( [T2] i, assuming Mx = 2 and re- 
placing n„ by using Eq. (17i. (iii) From the recombination 
order 



l+l = 



Eu 
kT 



1, 



(19) 



with an experimentally determined recombination order of 
2.6.^^ From (i) Kirchartz etal. found Eu w 50 meV, from 
(ii) £■[/ « 100 meV at room temperature, whereas from (iii) 
Eu =41 meV can be calculated. This discrepancy was dis- 
cussed in some detail in their publication.^ 

We considered a more general recombination rate, includ- 
ing the annihilation of free charge carriers with one another 
and free with trapped charge carriers. 



R{n) — k'nc{nc + «r) = k'rtcn. 



(20) 



Here, the overall charge carrier concentration is given as 
n = tic +nt- For simplicity, we used the same recombination 
prefactor k' for both contributions and 

Usually the rate limiting step in nongeminate recombina- 
tion is the finding of the actually localized charge carriers by 
a hopping process, which is reflected by a prefactor propor- 
tional to the mobihtyof the mobile charge carriers, based on 
Langevin theoryj ^ l ^ l I However, although beyond the scope 
of this article, we point out that in the multiple-trapping-and- 
release approximation of hopping transport, the effective mo- 
bility /J is proportional to a trap-free mobility /jq times the frac- 
tion of free to all charge carriers, /jq ■ Within this approach, 
R o= /urP' and R °^ /Joricti are in principle equivalent. We ask 
the reader to bear in mind that this representation can serve 
as a first approximation, but neglects the field dependence in 



a hopping model. For our purpose here, however, Eq. ( 20 1 
is sufficient to describe the recombination process in terms of 
the ideality factors discussed above. 

The small-signal method TPV yields the effective lifetime 
of all charge carriers, T„. For charge extraction, we as- 
sumed that all charge carriers were extracted, in accordance 
with Kirchartz et alP^ Thus, our experimentally determined 
effective recombination rate is R = n/z„ = {tic +«f)/T„ (c.f. 
Eq. (fTO])). Comparing our assumption for the recombination 



rate, Eq. ( 20 1 to our effective recombination rate, we find 



T„ = (k'tlc) 



-1 



(21) 



In order to test our assumption, we considered the voltage 
dependence of the constituents « and T„ of the effective re- 
combination rate by investigating their ideality factors ti„ and 
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rix, respectively. Using the exponential tail modell^ we also 
followed the three different routes outlined above. Calculat- 
ing the characteristic tail energy (i) directly from n„, we deter- 
mined Eu = 45 meV at room temperature. By route (ii), using 



our Eq. ( 12 1 instead of Eq. (18 1, we again found 45 meV with 
either riR or (illuminated) riid. In contrast to Kirchartz etal.^, 
our result from (ii) is consistent with (i). We point out that 
the authors had a rather high ideality factor of 1 .6 which was 
determined under dark conditions, in contrast to our consis- 
tent values of n,v/ ^ hr — 1.20. When using our dark ideality 
factor of 1.30, we find 48 meV for the characteristic tail en- 
ergy. Deriving the tail energy from the recombination order, 
route (iii). 



X+l = ^ + l 

tlx 

2__ 
kT rii 



(22) 
(23) 



using Eq. ( 17 i, yielded again a result compatible with the 



other two approaches, as essentially always the same equa- 
tions are used for all routes within our consistent framework. 



In contrast to Eq. ( 18 i, as used in Ref. 5, deviations of from 
the value of 2 are considered. We point out that due to our 
focus on the interpretation of ideality factors here, which is 
only qualitative in view of the recombination rates, we are not 
able to make statements concerning the contributions from de- 
layed recombination.'^^** For understanding the consistency 
of the charge carrier concentration and voltage dependence of 
the loss current in view of recombination order and ideality 
factor, see Appendix]] 

The increase of n„ with lower temperatures (Fig.|5](a)) also 
indicates a thermally activated process concerning the charge 
carrier concentration: the lack of thermal energy leading to a 
growing fraction of trapped charge carriers is compatible with 
the assumption of exponential tail states. Calculating the char- 
acteristic Urbach energy from the temperature dependent n„ 
by using Eq. ( 17 1, we found Eu — 40 meV for 200 K, show- 
ing that an exponential tail may not be the precise shape of the 
density of trap states, but can serve as an approximation in a 
limited temperature range. 



Within our assumption of the recombination rate, Eq. ( 20 1, 
the effective lifetime is inversely proportional to the free 
charge carrier concentration tic (Eq. (21 1). Therefore, we 
would expect that n-^ should equal 2 in accordance with «„ 
for free charge carriers. However, this statement is only valid 
if the recombination prefactor k' is assumed to be voltage in- 
dependent. While we cannot determine the detailed reason 
for our experimental finding of = 1.84, we point out that 
the assumption of k' ^ k'iy) may not hold tru d^^ l ^^ l and that 
concentration gradients are disregarded. 

In order to verify that the nongeminate recombination 
mechanism in annealed P3HT:PCBM solar cells is indeed due 
to losses of free carriers with one another and with carriers 
trapped in exponential tail states with Eu « 50 meV, we recon- 
sidered thermally stimulated current measurements presented 
previously.^" In this set of experiments, we did not account for 
trap states deeper than 400 meV. The distribution of the trap 
states is shown in a semilogarithmic plot (Fig.|6]l. Neglecting 
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FIG. 6. (Color Online) Trap density distribution of a P3HT:PCBM 
blend device, as measured by the thermally stimulated currents 
technique.-"^ Disregarding the detailed shape of the trap distribution, 
we approximated the energetic tail (solid line), finding a characteris- 
tic energy Eu ~ 57 meV. 



the detailed shape, we approximated the energetic tail by an 
exponential Urbach fit. Our analysis yielded a characteristic 
energy Eu ~ 57 meV, in accordance with a recent reconstruc- 
tion of the density of trap states in P3HT:PCBM blends,^' and 
verifies our findings from above. 



V. CONCLUSION 

The SE parameters such as ideality factor n,v/ and dark 
saturation current jo were derived either from static j/V- 
measurements or from the transient techniques (TPV/TPC), 
were shown to coincide within experimental error for OSC 
based on P3HT:PCBM. We directly compared the ideality fac- 
tor and dark saturation current density experimentally and re- 
produced the saturation photocurrent for P3HT:PCBM OSC 
over the temperature range from 200 to 300 K at various light 
intensities. Recombination current determined under open cir- 
cuit conditions is shown to be equal to the respective satu- 
ration current, which implies that the polaron pair dissocia- 
tion is not significantly influenced by the electric field. A 
good agreement of static and transient approaches over a wide 
temperature range demonstrates the validity of the Shockley 
model for OSC, if the charge carrier photogeneration is volt- 
age independent. Additionally, the dark saturation current 
measured at different temperatures was used to determine the 
effective band gap of P3HT:PCBM blend to be in the range 
of Eg K, 0.9 — 1.1 eV, which is in good agreement to values 
from literature. Considering the ideality factors, we found that 
nongeminate recombination of free with both free and trapped 
charge carriers in tail states is the dominant loss mechanism. 
Using data from thermally stimulated current measurements, 
we verified that the charge carrier traps can indeed be approx- 
imated by an exponential trap distribution. 
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Appendix: Voltage and charge carrier concentration 
dependence of the loss current 



In Eq. (16) of Ref. 5 Kirchartz etal. considered the carrier 
concentration dependent loss current at Voc, 
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and its connection to the voltage dependence of the carrier 
concentration, c.f. our Eq. (|4]). 
They found that 



din j I ass d\nn d In jioss 



d\nn 



dV 



kT 



dV 

"id kT 



(A.2) 



has to hold for self consistency (in our notation). Evaluating 



this equation, considering A, = ^, we see that 



(A.3) 



directly yields Eq. (12 1 
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